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Abstract 

FeGas is an unusual intermetallic semiconductor that presents intriguing magnetic responses to 
the tnning of its electronic properties. When doped with Ge, the system evolves from diamagnetic 
to paramagnetic to ferromagnetic ground states that are not well understood. In this work, we 
have performed a joint theoretical and experimental study of FeGas-ajGea; using Density Functional 
Theory and magnetic susceptibility measurements. For low Ge concentrations we observe the 
formation of localized moments on some Fe atoms and, as the dopant concentration increases, a 
more delocalized magnetic behavior emerges. The magnetic configuration strongly depends on the 
dopant distribution, leading even to the appearance of antiferromagnetic interactions in certain 
configurations. 

PACS numbers: 71.27.+a, 75.50.Bb, 71.15.Mb, 76.80.+y 
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I. INTRODUCTION 


The physics of Fe-based intermetallic semiconductors such as FeSi, FeSb 2 and FeGaa can 
present very unusual phenomena, creating an important framework for the development of 
condensed matter physics. These materials can present properties of strongly correlated 
electron systems, owing to the presence of the Fe d levels, concurrent with semiconducting 
gaps due to the strong hybridization between the Fe d and the post-transition metal or 
metalloid s,p levels. They have been studied mainly for their thermoelectric propertie^^. 
However, they are also interesting from fundamental physical aspectP', still presenting im¬ 
portant unsolved puzzles. 

Here we address the puzzling magnetic response upon electronic doping of FeGas, one 
of the few Fe-based materials known to be diamagnetic. This diamagnetism of the pure 
compound has also been assigned to the strong hybridization between the Fe 3d levels 
and the Ga 4s and 4p levels. Starting from this already unusual state, strikingly different 
magnetic responses can be observed upon chemical substitutions that lead to doping with 
holes or electrond^. 

There is much controversy in the literature regarding the most stable magnetic configura¬ 
tions of pure and doped FeGas, as well as on the mechanisms involved. Some works propose 
an antiferromagnetic ground state for pure and doped FeGagP^, whereas others report a 
diamagnetic ground statd^®. There are also authors invoking an itinerant origin for the 
ferromagnetic ground statd^, whereas other reports show a localized nature for the magnetic 
moment^^^^, or both!^. 

In order to understand this system, we have performed a joint theoretical and experi¬ 
mental work using First-principles Density Functional Theory (DFT) calculations and bulk 
magnetization techniques on Ge-doped FeGas single crystals. Our results confirm that, upon 
doping, the system quickly evolves to a mainly ferromagnetic (FM) ground state. However, 
the evolution of magnetic response to Ge doping is far from trivial. The magnetic moments 
are not evenly distributed throughout all Fe atoms as expected in an itinerant picture, and 
different groupings of magnetic moments on the Fe atoms are resolved, that depend on both 
doping concentration and dopant distribution. Additionally, some distributions do indeed 
lead to Fe moments aligned antiparallel to the mostly FM state. Our results point to an evo¬ 
lution from localized to delocalized character and an intricate competition between FM vs. 
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antiferromagnetic (AFM) pictnres for the magnetic response in the doped materials, which 
may lead to the conclnsion that almost all of the previons works were providing partially 
correct but incomplete descriptions of this fascinating system. 

II. METHODOLOGY 

Crystal growth and experimental characterizations 

Polyhedral single crystals of FeGa 3 _a;Ge 3 ; with x = 0.27 were grown by a standard Ga 
self-flux methocP^*^. High purity elements were sealed inside evacuated silica ampoules 
with a Fe:Ga:Ge proportion of l:[3(3-x)]:[3x]. These ampoules were heated to 1100 °G and 
cooled to 550 °G over 150 h, then removed from the furnace and quickly centrifuged for 
separation of the molten flux. The effective doping level of each crystal was estimated by 
comparing the Sommerfeld coefficient, transition temperature and effective moment with the 
results by Umeo et al^. X-ray diffraction patterns of powdered crystals were obtained on a 
Bruker D8 Focus machine and are all consistent with the FeGaa type structure, with rehned 
lattice parameters of the pure FeGas in agreement with published work^^^Il jyQ magnetic 
susceptibility for temperatures between 2 and 300 K was measured on a Quantum Design 
SQUID-VSM. For applied helds below 30 Oe the flux gate accessory was used to cancel any 
remnant held in the superconducting magnet, then a sequence of zero held cooled (ZFG), 
held-cooled cooling (FGG) and held cooled warming (FGW) branches was performed. 

Density Functional Theory calculation 

Spin-polarized hrst-principles calculations based on DFT were carried out using the 
Generalized Gradient Approximation with a parametrization targeted especially for solids, 
namely PBEsoP^. Also, we performed calculations introducing an on-site Goulomb interac¬ 
tion term, the DFT-I-U method^^*^, characterized by the Hubbard U parameter acting on 
Fe d-states. The Kohn-Sham equations were solved using the Projected Augmented Wave 
(PAW) method^, as implemented in the Vienna ab-initio Simulation Package (VASP)^^^. 
The PAW atomic reference conhgurations are 3p®3d®4s^, 3d^°4s^4p, 3d^°4s^4p^ for the Fe, 
Ga and Ge respectively. These settings reproduce correctly the properties of this material 
as reported in Ref. HI The kinetic energy cut-oh for the plane-wave expansion was 400 eV. 
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For each dopant configuration, the lattice parameters were fixed at the FeGaa equilibrium 
crystal structure while the atomic positions were allowed to relax without any symmetry 
constraint. In the relaxation process, Monkhorst-Pack k-point meshes were used, and de¬ 
pending of the size of the supercell we used F-centered k-point meshes of 6 x 6 x 6, 6 x 6 x 3, 
and 2x2x2 for the calculation of charge and magnetization densities. 

III. RESULTS AND DISCUSSION 

FeGaa has a tetragonal crystal structure belonging to the space group PA 2 /mnm (136) 
with one Fe atom located at the Wyckoff position 4 /(m, m, 0) and two types of Ga atoms: 
Gal occupies the Wyckoff position 4c(0, ^,0) and Ga2 is located at the Wyckoff position 

The simplest way to theoretically model n-type doping in FeGas is by shifting the Fermi 
level of the system up, which will represent the occurrence of electrons in the conduction 
band. This preliminary approach leads to the observation that the magnetic response de¬ 
pends on the electron concentration on the system. For low electron densities (< 0.09e/f.u.) 
we do not observe any perturbation of the system, and no magnetic moments are observed 
in the Fe atoms. As we increase the electron density, we start to observe a magnetic moment 
in some Fe atoms, i.e., the electrons are not homogeneously distributed through the whole 
lattice. In this case, the electrons seem to be self-trapped in a lattice distortion. A similar 
self-trapping behavior has been reported before for the semiconductor MnCP^l. At higher 
electron concentration (0.25e/f.u.) it is observed a more delocalized (itinerant) magnetic 
state with similar magnetic moments in all Fe atoms (see Fig. [^. Test calculations have 
shown that these results are robust with respect supercell sizes. 

After modeling the system with a preliminary effective doping model, we advance our 
calculations to an explicit doping approach of FeGas with Ge atoms (FeGa 3 _a:Gea:). For 
this study, the host system is seen as a layered structure with the Gal ions placed at the 
same plane of Fe and Ga2 occupying its own separate plane. Galculations were performed 
at concentrations x = 0.03, 0.06, 0.09, 0.13, 0.19, 0.25, 0.38 and 0.50 for several different 
distributions of substitutional Ge impurities throughout the supercell, always on Ga sites, 
with the appropriate amounts to obtain the desired concentrations. We have also converged 
different spin configurations for each atomic distribution, including (i) non-spin polarized. 
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(ii) all Fe atoms with parallel spin and (iii) all Fe dimers with antiparallel spins. Along this 
work we will discuss the most stable spin conhguration at each case. 

When a single Ge atom is inserted in a 2 x 2 x 2 supercell at either Gal or Ga2 site 
(x = 0.03), the Ga2 site is preferred by an energy difference of 0.11 eV. In this case, the Ge 
atom also induces a tiny magnetic moment on its Fe neighbors (5 x 10“^ps/f.u.), whereas 
at the Gal site it does not lead to any magnetic moment on its Fe neighbors. This contrast 
demonstrates the complex nature of this material, since the exact lattice position where 
each Ge atom enters directly influences the magnetic conhguration, and actual samples will 
inevitably feature a statistical distribution of the dopant among the nearest neighbors of an 
Fe atom. Gonsequently, a complete description of the magnetic behavior can only be achieved 
through a meticulous and thorough exploration involving many different distributions, which 
we undertake next. 

By increasing the dopant concentration in different supercells, we were able to obtain 
several different magnetic conhgurations, wherein the total moment per f.u. also resulted 
different. For a given concentration, there can be different spin conhgurations that diher 
simply by the lattice distribution of the Ge atoms. This is exemplihed in the main panel 
of Fig. [^a), where the total magnetic moment per f.u. is plotted for several diherent 
distributions of a hxed dopant concentration. The blue diamonds in this hgure indicate 
the most stable spin conhguration. For x = 0.03 and 0.06 there exist distributions with 
null magnetic moment. For each concentration and distribution, we also obtain diherent 
values of magnetic moment on the Fe atoms, depending on their positions with respect to 
the dopants. This result is analyzed by counting the number of Fe atoms with a certain 
magnetic moment in each conhguration, then applying a Gaussian distribution around each 
value as a simplihed representation of conhguration distributions in actual samples. The 
inset of Fig. [^a) shows the resulting magnetic moment distribution on the Fe atoms for 
the lowest energy conhgurations with x = 0.09, 0.38 and 0.50. For x = 0.09, the largest 
peak (generated by 22 Fe atoms) is around zero magnetic moment. Five Fe moments take 
values around 0.15/1^ and the others take values around -0.17, -0.12, 0.22, 0.30 and 0.33 /is. 
Note that the lower concentrations have a wider distribution of peaks, whereas the higher 
concentration shows a clear trend towards narrowing the distribution around a single value. 
These results clearly pose a question on the validity of any calculation performed within the 
virtual crystal approximation, since through such a method only an average ehect can be 
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probed and all the local site information is lost. 

We have also performed calcnlations adding the Hnbbard U term to the energy fnnctional 
(DFT+U) within the simplihed rotationally invariant approachP^ using [4// = 2 and 3 eV 
for the 2 X 2 X 2-supercell. Figure [^b) shows the magnetic moment distribution on the Fe 
atoms when one Ge impurity occupies either Gal or Ga2 positions. The inclusion of the 
Hubbard term induces a local moment on some Fe atoms that is quantitatively incorrect in 
comparison to the experimental results. This proves that the use of a semi-local functional is 
the correct approach to calculate the physical properties of this weakly correlated material. 

With the aim of visualizing the real-space distributions of magnetic moments in the 
system, we plot the calculated spin densities for the more stable distribution of Ge atoms at 
each concentration (Fig. |^a)-(h)). In these representations one can observe that the induced 
magnetic moments (yellow surfaces representing spin densities in an arbitrary 2 ;-direction) 
are initially localized close to the impurity as we can see in Figs, [^a)-(d). By increasing 
the dopant concentration x, the magnetization gradually spreads throughout the Fe atoms 
on the lattice, leading to a more uniform and delocalized character of the magnetization 
at high concentrations as shown in Figs. |^f)-(h). Also, our results show that the induced 
magnetization does not depend on the supercell size, this is clearly observe in Fig. [^f) (x = 
0.25) where the magnetization densities for the unit cell, 1x1x2 and 2x2x2 supercells 
are fairly similar. It is also worthy of notice that some of the lowest energy conhgurations 
exhibit the formation of a few antiparallel spins on the Fe dimers, exemplified in Fig. |^c) 
and (e) by the blue surfaces seen for x = 0.09 and 0.19. In most of the distributions the 
induced magnetic moments are oriented in the same direction, leading to fully ferromagnetic 
states. Thus, two main points can be derived from the PBEsol calculations which led to 
the results in Fig|^ (i) an initial emergence of localized moments on some Fe atoms nearby 
diluted Ge impurities evolves, as a function of dopant concentration, into a more delocalized 
character at larger concentrations, ending in similar moments on all Fe atoms; (ii) in certain 
conhgurations it becomes energetically favorable for a few Fe atoms to align antiparallel to 
their neighbors, indicating the presence of a minor AFM component over the majority FM 
interaction. 

For experimental support of the above simulation results, we performed magnetic suscep¬ 
tibility measurements on a Ge-doped single crystal with x = 0.27, which places it far enough 
into the magnetically ordered ground state region to avoid issues related to non-Fermi-liquid 
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behavioil^. Fig. |^a) presents the T-dependence of the DC magnetic susceptibility of the 
FeGa 2 . 73 Geo .27 crystal under external magnetic held of FT = 10 Oe for the ZFG (black), FGG 
(red), and FGW (blue) branches. In the ZFG branch the response in the saturated region 
(0.30 emu/mol) is lower than the other branches, and shows a small initial increase up to 
about 20 K. This may be the result of the formation of magnetic domains in the single 
crystal upon cooling at zero held, which are not realigned at 2 K due to the application 
of such a weak magnetic held. Increasing temperatures then allow some of the domains 
to realign themselves. In the subsequent FGG branch, there is a small but well dehned, 
anomalous peak in the magnetic susceptibility around 48 K in which the susceptibility in¬ 
creases sharply upon cooling through Tc, then decreases equally sharply a little but returns 
to a slowly increasing behavior upon further cooling below 40 K, eventually reaching 0.33 
emu/mol 2 K. In this case both the magnetic domains formed at Tc and the magnetic mo¬ 
ments were already under an applied held which dehned a clear symmetry breaking direction 
for parallel and anti-parallel alignment. Very similar behavior is visible in previous works 
reporting only FGG branched but was not discussed. This rare behavior disappears in the 
subsequent FGW branch, that gradually recovers similar behavior to the ZFG branch as the 
sample is warmed through Tc- Surprisingly, a similar behavior appears in a weak ferromag- 
net compound MnSi classihed as a helimagnelP^. This interesting feature is expected for 
non-centrosymmetric structures, which is not our case. 

In order to further explore the behavior observed in the FGG curve, measurements under 
gradually higher applied magnetic helds {H = 30 Oe, 100 Oe, 300 Oe, FGG mode) were 
carried out and are shown in Fig. |^b). The peak is still clearly visible at if = 30 Oe and 
detailed in the inset of Fig. |^b). Under 100 Oe the FGG peak has practically disappeared, 
and under 300 Oe the sample recovers the expected shape of FGG curves for traditional 
ferromagnets, i.e., a sharp increase around Tc followed by further, gradual increase towards 
saturation as thermal energy is removed from the spin system. The small but sharp decrease 
just below Tc observed in low-held FGG curves is therefore interpreted as evidence of a 
minority of antiferromagnetic interactions present in the sample, resulting in some of the 
Fe moments aligning antiparallel to the majority, as was found in the DFT calculations for 
X = 0.09 and 0.19 lowest energy conhgurations (Figs. I^c) and i e)). Although we cannot 
rule out the inhuence of magnetic domains, we argue that it is only possible to see this weak 
antiferromagnetic contribution when all of the domains and spins organize under a well- 
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defined alignment axis (FCC case). Also, this is clearly not a trivial case of ferrimagnetism, 
since the antiparallel moments align themselves at a discernibly lower temperature than the 
majority (possibly indicating a certain level of frustration) and because the antiferromagnetic 
component proves to be quite fragile. For higher external fields the peak quickly disappears 
and the susceptibility curves adopt the more traditional, monotonic increase upon cooling 
that is expected of a ferromagnet. It is worth noting that the different values of susceptibility 
saturation at low temperature for these three applied magnetic fields are still compatible 
with an almost linear behavior of the magnetization curve at low fields, within the resolution 
of our experiments. By fitting the inverse susceptibility with a modified Curie-Weiss law at 
high temperatures (see Fig. |^, we extract /ig// = 0.84(6)/rB/f-u, 9cw = 56(3) K and xo = 
3.7(8) X 10“^emu/mol. 

The inset of Fig. [^presents a magnetization measurement showing an almost linear be¬ 
havior at low fields and no significant hysteresis within the precision of the SQUID measure¬ 
ment. In the fully ferromagnetic state, the saturated magnetic moment ^sat = 0.19(l)/iB/f.u. 
at T = 2 K is comparable to the theoretically calculated average of 0.22 ps/Fe for the more 
stable configuration with x = 0.25, (see Fig. [^a)). It is worth noting that the effective 
moment, fieff, is different compared with the saturation moment, fisat- The itinerant vs. 
localized character of ferromagnets can be characterized by the Rhodes-Wolhfarth Ratio, 
RWR = For a localized system, the value of RWR should be close to 1, while 

it diverges for itinerant ferromagnets. In our case, RWR is about 4.4, supporting the sce¬ 
nario of intermediate behavior between localized and itinerant forms of the ferromagnetism 
as our DFT results showed. Finally, the evolution from localized towards more delocal¬ 
ized magnetic character as a function of dopant concentration, and also the coexistence of 
FM and AFM components, were also supported by recent Nuclear Quadrupolar Resonance 
measurement^^. 


IV. CONCLUSIONS 

In conclusion, we have shown that the magnetism of Ge-doped FeGas is in fact much 
richer than what had been described so far. Magnetic moment distributions throughout the 
Fe atoms are complex and strongly dependent on Ge dopant concentration, as well as on 
the manner in which they are distributed in the lattice. For small dopant concentrations. 


Fe site symmetry is broken and different types of Fe atoms emerge: some bearing localized 
moment (induced by nearby Ge) and the rest remaining non spin polarized. This scenario 
gradually evolves with increasing Ge concentration to a more delocalized character until 
the dominance of itinerant magnetism is established at the highest experimentally achieved 
doping levels {x ~ 0.41). Additionally, the magnetic ordering is found to be mainly FM, 
but in some cases it is found a minor AFM component brought in by Fe moments that align 
antiparallel under circumstances not yet fully understood (once again strongly dependent 
on the dopant concentration and configuration). These predictions made by our meticulous 
DFT calculations were very well supported by magnetic susceptibility measurements on a 
FeGa 2 , 73 Geo .27 crystal, that showed a fragile AFM component distinguishable from the major 
FM transition. 

What makes this fascinating system even more compelling is that all these rich features 
arise even before mixing in the broad range of anomalies that accompanies its reported quan¬ 
tum critical behavior around x = 0.13, which is not contemplated in either the theoretical 
modeling or the experiments, intentionally performed here on a sample that is well into the 
magnetically ordered region. 

Note added to arXiv version: this paper was originally submitted as a collaboration 
which included preliminary ®^Fe Mossbauer data, but after further consideration the authors 
involved in those experiments have chosen to prepare a separate, comprehensive publication 
to be submitted briefly. 
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(a) 0.06e/f.u. 



(b) 0.13e/f.u. (c) 0.25e/f.u. 


FIG. 1: Electron doped FeGaa magnetization density of a 2 x 2 x 2-supercell: (a) small 
electron density (0.06e/f.u.), (b) medium electron density (0.13e/f.u.), and (c) large 
electron density (0.25e/f.u.). Brown and grey spheres represents Fe and Ga atoms, 
respectively. The yellow surface represents the magnetization density around each atom. 


12 























(a) PBEsol 



FIG. 2: (a) PBEsol FeGa 3 _a;Gea; average magnetization per f.u. as a function of the Ge 
concentration x for different substitutional distributions. The blue diamonds indicate the 
lowest energy configuration for each concentration. Inset: Fe magnetic moment 
distribution of the lowest energy configurations for x = 0.09, 0.38 and 0.50. (b) DFT+U Fe 
magnetic moment distribution for a 2 x 2 x 2-supercell where one substitutional impurity 
is located in either Gal or Ga2 site (x = 0.03, Ueg = 2 and 3 eV). 
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(a) x=0.03 


(b) x=0.06 



(c) x=0.09 (d) x=0.12 (e) x=0.19 



(f) x=0.25 



(g) x=0.37 (h) x=0.50 


FIG. 3: FeGas-j^Gej; magnetization density for the lowest energy configuration of some Ge 
concentrations x and supercell sizes. Brown, gray and green spheres represents Fe, Ga and 
Ge atoms, respectively. The yellow and blue surfaces represent positive and negative 

values of the magnetization density. 
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FIG. 4: a) T-dependence of the magnetic susceptibility at if = 10 Oe in the ZFC, FCC 
and FCW modes for FeGa 2 . 73 Geo. 27 - b) T-dependence of the magnetic susceptibility at 
if = 30 Oe, 100 Oe and 300 Oe. The inset zooms in around the peak region that evidences 

an anomalous antiferromagnetic behavior. 
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FIG. 5: Inverse susceptibility versus temperature for FeGa 2 , 73 Geo. 27 - The red line is the 
fitting in high temperatures of the Gurie-Weiss law. The inset shows the magnetization 
(M) as a function of magnetic field (H) in T = 2 K. 
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